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Abstract 

We study CP violation and final state interactions between pions and kaons in B + , B~ , — B° and 
B° decays into Kirir. The weak transition amplitudes consist of two terms: the first part is derived 
in QCD factorization approach and the second one is a phenomenological long-distance charming 
penguin contribution. The final state Ktt interactions in S- and P-waves are described by strange 
scalar and vector form factors, respectively. These are determined using a unitary coupled channel 
model together with chiral symmetry and asymptotic QCD constraints. The final state interactions 
are dominated by presence of the scalar ^(1430) and the vector if* (892) resonances. We show 
that additional charming penguin amplitudes are needed to reproduce the latest experimental Ktt 
effective mass and helicity angle distributions, branching fractions and asymmetries obtained by 
Belle and BaBar collaborations. 

Introduction 

Recent detailed Dalitz plot analyzes for different B — > Ktt + tt~ decays were performed by BaBar 
and Belle collaborations [H [21 O [U [5] . New and rich data on Kir effective mass (niKn) and helicity 
angle distributions, branching fractions and CP asymmetries have been delivered. The effective mass 
distributions show distinct accumulation of events below 2 GeV around the scalar Kq(1A30) and 
the vector if* (892) resonances. In experimental analyzes the event distributions have been studied 
using the isobar model with decay amplitudes parameterized by sums of Breit-Wigner terms and a 
background. In [3l 0] an effective range component has been introduced in the Ktt 5-wave amplitude. 

In this short report we present main results of our work on B decays into three mesons, B —> 
Ktt + tt~, with final state interactions between kaons and pions in the S- and P-waves [6]. In our 
analysis we restrict ourselves to the low effective Ktt mass region rriKn < 2 GeV where most of the 
Ktt resonant structures are seen. It is reasonable to expect that in such a kinematical configuration 
three-body interactions are suppressed. This allows us to grant the validity of quasi two-body factor- 
ization approach. We also assume that the Ktt pairs originate from a quark-antiquark state. Main 
assumption in our approach is, however, that the corrections to factorization can be absorbed into 
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niK-n independent modifications of the QCD Wilson coefficients. We take parts of these corrections 
from quasi-two-body calculations and append a phenomenological part. 

In the factorization approach the decay amplitudes can be expressed as a product of the effective 
QCD part and the two matrix elements of the vector currents which involve two functions of m 2 KlT : the 
strange scalar and the vector form factors. The introduction of form factors, constrained by theory 
and experiments other than B decays, is an alternative to the isobar model often used in experimental 
analyzes. This latter approximation violates unitarity what can, for example, lead to a distortion of 
information on resonances present in the final states. 



2 Theoretical model 

The weak amplitudes for the B^, B° and B° decays, have one term derived in the QCD factorization, 
and the second one being a phenomenological contribution to the penguin amplitudes. For example, 
the Kit S-wave contribution to the B~ — > K~tt + tt~ decay amplitude reads (expressions for other B 
decay amplitudes can be found in [6]): 
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where the a, (i = 4,6,8,10) are the leading order factorization coefficients with 0(a s ) vertex and 
penguin corrections added in the quasi-two-body approximation of pseudoscalar-scalar or pseudoscalar- 
vector final states. The Cj coefficients are additional phenomenological complex parameters fitted to the 
experimental data. They can be interpreted as contributions from charming penguins, hard spectator 
interaction in the perturbative regime and as annihilation terms. Gf denotes the Fermi coupling 
constant, Mb, ttlk and are the masses of the charged B mesons, kaons and pions, /q n (m 2 Kn ) 
and n {m^w) are the &~ to 7r~ and the Ktt scalar form factors respectively and m&, m s and 
are b, s and d quark masses. The \ u = V u f,V* s and A c = V c bV* s are products of the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix elements V qq i. As we have shown in [B], the ma 2 KT . dependence of the 
S- wave amplitude (3rd line of the Eq. (1)) plays an important role in the description of the sizable 
enhancement in the Ktt mass distributions below 1 GeV. In the isobar model the corresponding m^n 
dependence is approximated by one fitted constant parameter. Such approximation may not be very 
good for the wide -K"g(1430) resonance since the value of m 2 Klx varies by a factor almost equal to 10 
from the Ktt threshold to the limit of our calculations about 1.8 GeV. 

The scalar and vector Kit form factors are connected to the scattering amplitudes in the S- and 
P-waves via unitarity relations. In the S-wave we treat two coupled Ktt and Krf channels and in the 
P-wave three ones: Ktt, K*tt and Kp. The corresponding scattering amplitudes are constrained by 
experimental data, especially by the LASS results obtained at SLAC, and by theoretical constraints 
at low effective mass. The Muskhelishvili-Omnes equations are used to calculate the Ktt form factors. 



3 Numerical results 

The Ktt effective mass and helicity angle distributions, branching fractions, CP asymmetries and the 
phase difference between the B° and B° decay amplitudes to K*(892)tt are calculated from our model 
and fitted, using the minimization program with four phenomenological free parameters c^'g, to 319 
data obtained by BaBar and Belle collaborations. Branching fractions play a special role as they 
determine the absolute size of the decay amplitudes. As an input we use the branching fractions 
for the B — > K*(892)tt decays which are well determined since the Ktt P-wave amplitude below 1 
GeV is entirely dominated by the K*(892) resonance. Contrary to the P-wave amplitude, the S- wave 



amplitude contains not only the resonant part related to the wide -fTg(1430) but also a very broad 
i^g(800) below 1 GeV. This complexity of the S- wave results in different parameterizations of the 
amplitudes in the BaBar and Belle analyzes. Therefore we do not use the experimental branching 
fractions for the B — > i£g(1430)7r, which are not well determined, being largely model dependent. 
Some of our model predictions for branching fractions in both waves are given in Table 1. These 
numbers have been calculated by integration of the double differential branching fractions in limited 
ranges of vtik-k around the -K"q(1430) and X*(892) resonances. The model branching fractions for the 
B —7- if*(892)7r decays agree quite well with the corresponding experimental values within their errors. 
If the phenomenological parameters cf (i = 4, 6 and p = u,c) are equal to zero, than the theoretical 
branching fractions for the P-wave are 4-5 times smaller than the experimental ones and those for the 
S-wave differ by a factor of 2. One can notice that the theoretical results for the S"-wave are lower 
than the experimental data, being substantially smaller than the Belle results but closer to the BaBar 
numbers. 

Table 1: Branching fractions averaged over charge conjugate B — > Kirn decays in units of 1CP 6 . 
(K Jr -K~)p and (K + ir~)s denote the (K + ir~) pair in P-wave and S"-wave, respectively. The values of 
the model, calculated by the integration over the given uik-k range, are compared to the correspond- 
ing Belle [U [2] and BaBar [3j H] results. The model errors are the phenomenological parameter 
uncertainties found in the minimization procedure. 



decay mode 


rriK-K range (GeV) 


Belle 


BaBar 


our model 


B + (K+tt-)p ir- 


(0.82,0.97) 


5.35 ±0.59 


5.98 ± 0.75 


5.73 ±0.14 


B° -> (^ U vr+) P 7T" 


(0.82,0.97) 


4.65 ±0.77 


6.47 ± 0.75 


5.42 ±0.16 


B + (K + ir-) s tt- 


(0.64,1.76) 


27.0 ±2.5 


22.5 ±4.6 


16.5 ±0.8 


£ u -»• (^ u vr+) s 7T- 


(0.64,1.76) 


26.0 ±3.4 


17.3 ±4.6 


15.8 ±0.7 



In Table 2 we present comparison of our results for the CP asymmetries with experimental ones. 
For B~ (K~ir + )pir~ decays our asymmetries lie between those of Belle and BaBar. The results 
for B~ — > (K~it + )s 7r~ and for B° — > (K°7r~)s 7r + decays agree with the experimental values of both 
collaborations. 

Table 2: Direct CP asymmetries averaged over charge conjugate reactions. The values of the model, 
calculated over the indicated rriKn range, are compared to the Belle [1] and BaBar [3j 0] results. The 
model errors are the phenomenological parameter uncertainties found in the minimization procedure. 



Decay mode 


rriK-w range (GeV) 


Belle 


BaBar 


our model 


B~ — > (K~ir + )p7r~ 


(0.82,0.97) 


-14.9 ±6.8 


3.2 ±5.4 


-2.5 ± 1.3 


B {) -> (K U TT-) P 7T+ 


(0.82,0.97) 




14 ±12 


-19.6 ±3.0 


B~ (K-ir+)sn- 


(1.0,1.76) 


7.6 ±4.6 


3.2 ±4.6 


5.4 ± 1.0 


B u -> (K u vr-) s vr + 


(1.0,1.76) 




17 ±26 


-0.2 ±1.3 



The fit on the model parameters leads to a good agreement with the experimental data, particu- 
larly for the kaon-pion effective mass and helicity angle distributions. Examples of these results are 
presented in Figs Q] and [2j 

In Fig. [3]we show decomposition of the averaged mjcn distributions of the charged B meson decays 
into a resonant, background and interference terms. We compare these components calculated using 
our amplitudes and the BaBar's LASS parameterization. The background part in our approach is 
several times smaller than that effective range term in experimental analyzes. Worthy of note is the 
fact that interference parts have very different sizes and even different signs after integration over 
rriKir- 
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Figure 1: The J fT ± 7r T effective mass distributions in the B — > K^tt^tt^ decays from the fit to the 
experimental data pp. The dashed line corresponds to the S- wave contribution of our model, the 
dotted line that of the P-wave and the histogram represents the coherent sum of the S- and P-wave 
contributions. 

4 Conclusions 

We have shortly presented main results of our analysis of the Kir interactions in the B — > Ktt + tt~ 
decays with Ktt pairs in the S- and P- waves. In our approach sufficiently good description of ex- 
perimental data can be achieved after introduction of four complex phenomenological coefficients into 
amplitudes. These additional parameters are added to the QCD factorization coefficients and can 
be interpreted as coming from long-distance charming penguin amplitudes, hard spectator and weak 
annihilation contributions. 

We have shown that the Ktt P-wave amplitude is entirely dominated by If* (892) whereas the S- 
wave with the wide resonance Kq (1430) and the very wide enhancement below 1 GeV (often interpreted 
as the -?Tq(800) meson) is more complex. In order to correctly account for these resonances we use 
scalar and vector form factors in the S- and P-waves instead of the Breit-Wigner amplitudes, often 
employed in the isobar model. This leads to the background term in the effective mass distributions 
which is several times smaller than the one obtained in each experimental analysis. 

We therefore propose to use in the future experimental analyzes the following parameterization of 
the S- wave B — > (Ktt)stt decay amplitude: 

A s =f^(m 2 K7r )(c /m 2 K7r + c 1 ), (2) 

where f^™ is the scalar Ktt form factor, while cq and c\ are complex numbers to be fitted from the 
data. Numerical values of the complex scalar form factor can be provided on request. 

Let us emphasize at the end that the amplitudes derived in our analysis [6j go beyond the two- 
body approach applied to such decays as B — > K*ir. We have proved that they correctly incorporate 
the two body final state interactions between K and ir if the effective Kn mass in the three-body 
B — > Ktt + tt~ decays is limited to about 2 GeV. 
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Figure 2: Helicity angle distributions for K^ir^ir^ decays calculated from the averaged double 

differential distribution integrated over ra^±^ mass from 0.82 to 0.97 GeV in the K* (892) case a) and 
from 1.0 to 1.76 GeV in the -K"q(1430) case b). Data points are from Ref. [5]. Dashed lines represent 
the S'-wave contribution of our model, dotted lines that of the P-wave and the dot-dashed lines that 
of the interference term. The histograms correspond to the sum of these three contributions. 
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Figure 3: Different components of the averaged rriKu distributions of the — > [K^t:^)s^ de- 
cays for: a) our model and b) BaBar's LASS parameterization [3]. The dashed lines correspond to 
the resonant _K"o(1430) contributions, the dotted-dashed lines to the background, dotted lines to the 
interference and the solid lines to their sum. 
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